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By  selecting  two  dressed  rotational  states  of  ultracold  polar  molecules  in  an  optical  lattice,  we  obtain  a 
highly  tunable  generalization  of  the  t-J  model,  which  we  refer  to  as  the  t-J-V-W  model.  In  addition  to  XXZ 
spin  exchange,  the  model  features  density-density  interactions  and  density-spin  interactions;  all  interactions 
are  dipolar.  We  show  that  full  control  of  all  interaction  parameters  in  both  magnitude  and  sign  can  be 
achieved  independently  of  each  other  and  of  the  tunneling.  As  a  first  step  towards  demonstrating  the 
potential  of  the  system,  we  apply  the  density  matrix  renormalization  group  method  to  obtain  the  ID  phase 
diagram  of  the  simplest  experimentally  realizable  case.  Specifically,  we  show  that  the  tunability  and  the 
long-range  nature  of  the  interactions  in  the  t-J-V-W  model  enable  enhanced  superfluidity.  Finally,  we  show 
that  Bloch  oscillations  in  a  tilted  lattice  can  be  used  to  probe  the  phase  diagram  experimentally. 
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Experiments  with  ultracold  atoms  have  recently  ex¬ 
tended  the  range  of  candidate  systems  for  realizing  uncon¬ 
ventional  states  of  matter  and  enabled  the  simulation  of 
models  describing  condensed  matter  phenomena  [1].  One 
major  goal  of  current  research  at  this  interface  between 
condensed  matter  and  atomic  physics  is  to  emulate  the 
Heisenberg  and  t-J  models,  which  are  believed  to  under¬ 
lie  certain  quantum  magnetic  materials  [2]  and  high- 
temperature  superconductors  [3],  respectively.  However, 
in  the  ultracold  atom  realization  of  these  models,  the  small 
superexchange  interaction  J  [1]  makes  the  underlying  phys¬ 
ics  extremely  challenging  to  observe.  At  the  same  time, 
ultracold  polar  molecules  [4],  such  as  KRb  [5,6]  and  LiCs 
[7],  have  recently  been  produced  in  their  electronic  and 
rovibrational  ground  states.  In  this  Letter,  we  show  that 
when  such  molecules  are  localized  in  an  optical  lattice,  their 
rotational  degree  of  freedom  can  be  used  to  simulate  tuna¬ 
ble  Heisenberg-like  models  at  unit  filling  of  the  lattice  and, 
in  the  presence  of  doping,  a  fully  tunable  generalization  of 
the  t-J  model  that  we  refer  to  as  the  anisotropic  t-J-V-W 
model.  Dipole-dipole  interactions  that  give  rise  to  this 
model  are  orders  of  magnitude  stronger  than  superexchange 
interactions  J  in  ultracold  atoms  and  can  therefore  better 
compete  with  other  relevant  energy  and  time  scales  such  as, 
for  example,  those  responsible  for  decoherence.  Moreover, 
we  show  that  the  resulting  long-range  interactions  are  fully 
controllable  with  dc  electric  and  continuous-wave  micro- 
wave  fields.  As  a  first  step  towards  demonstrating  the 
potential  of  the  model,  we  use  the  density  matrix  renormal¬ 
ization  group  method  [8]  to  obtain  the  ID  phase  diagram  for 
the  simplest  experimentally  relevant  case  and  show  that,  at 
low  fillings,  the  superfluid  phase  is  enhanced  relative  to  the 


one  in  the  conventional  t-J  model  [9].  We  propose  to 
probe  the  phase  diagram  by  using  center-of-mass  Bloch 
oscillations.  Given  that  KRb  has  already  been  loaded 
into  a  3D  lattice,  our  proposal  is  applicable  to  current 
experiments. 

The  Hamiltonian  and  its  features. — We  consider  dia¬ 
tomic  polar  molecules  in  their  electronic  and  vibrational 
ground  state  partially  polarized  by  a  dc  electric  field  along 
z,  confined  to  the  x-y  plane  [10],  and  loaded  in  that  plane 
into  the  lowest  band  of  a  square  optical  lattice.  As  de¬ 
scribed  below,  microwave  fields  are  used  to  isolate  in  each 
molecule  two  dressed  rotational  states  | m0)  and  \mf)  and  to 
obtain  the  t-J-V-W  Hamiltonian 


H  =  -t  S  i>«L<>  +  H-c-]  +  -  R,-i-3 

{i,j)m  14=  j 

r./,  j,  v  i 

X  L  2  SJS-  +  -^S)S)  +  -ninj  +  WntS)  J. 


(1) 


The  two  terms  describe  tunneling  and  dipole-dipole  inter¬ 
actions,  respectively;  ()  denotes  nearest-neighbor  bonds. 
Specifically,  cjm  creates  a  fermionic  [5]  or  bosonic  [6,7] 
molecule  on  site  j  (position  R  •)  in  dressed  rotor  state  m  E 
{m0,  m  | };  in  this  Letter,  we  focus  on  fermions.  Large 
reaction  rates  [10]  between  two  molecules  on  the  same 
site  enforce  the  hard-core  constraint.  We  define  njm  = 

nj  Sj  jm  \ ’  jm0 

njmi)/2  and  use  units  in  which  h  =  1.  The  Jz,  V,  and  W 
terms  can  be  understood  by  thinking  of  \m0 )  and  \mf)  as 
classical  permanent  dipoles  oriented  along  z,  while  the 
|  m0)  —  |  m  | )  transition  dipole  moment  gives  rise  to  the 
J y  term.  We  tune  these  dipole  moments  by  constructing 
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dressed  states  |  m0)  and  |  m , )  out  of  bare  rotor  states  using 
microwave  fields.  This,  in  turn,  allows  for  the  full  controll¬ 
ability  of  Jz,  J±,  V,  and  IV,  which  is  one  of  the  main  results 
of  the  present  Letter.  Lattice  Hamiltonians  based  on  more 
than  one  molecular  rotational  state  have  been  considered 
before  in  Refs.  [1 1-22].  An  important  difference  of  Eq.  (1) 
from  the  Hamiltonian  studied  in  Ref.  [17],  which  is 
most  closely  related  to  our  work,  is  the  presence  of  the 
Ji  term  [1 1], 

Equation  (1)  possesses  several  aspects  that  can  lead  to 
unconventional  many-body  phases.  Most  notably,  all  the 
interactions  are  long-range;  in  particular,  repulsive  long- 
range  density-density  interactions  are  expected  to  stabilize 
superfluid  correlations  [23].  Furthermore,  the  W  term  [24] 
can  break  SU(2)  symmetry  of  the  system  even  when 
Jz  =  J It  is  also  crucial  that  Jz  and  can  be  tuned  in 
both  sign  and  magnitude  (up  to  —100  kHz  in  LiCs)  inde¬ 
pendently  from  each  other  and  from  t.  This  contrasts  with 
the  cold  atom  realization  of  the  t-J  model,  where  7  «  /. 
Finally,  we  note  that,  in  the  limit  of  unit  filling  (n  =  1), 
only  the  terms  with  Jz  and  J ±  survive,  realizing  an  XXZ 
model  with  long-range  dipolar  interactions  and,  for 
Jz  =  7 x,  the  Heisenberg  model. 

Tuning  Jz,  J±,  V,  and  W. — We  now  show  how  to  use  dc 
electric  and  microwave  fields  [12,14-18,22,25-29]  to 
control  Jz,  J±,  ^  and  W,  leaving  a  detailed  derivation 
of  Eq.  (1)  to  Ref.  [30].  A  polar  molecule  in  a  dc  electric 
field  is  described  by  the  rigid-rotor  Hamiltonian  H{]  = 
B N2  —  d(]  E,  with  the  rotational  constant  B,  the  angular 
momentum  operator  N,  and  the  dipole  moment  operator 
d.  We  set  dp  =  ep  •  d,  where  e0  =  z  and  e±  = 
+  (x  ±  iy)/V2  [26].  At  E  =  0,  H0  has  eigenstates  | N,  M) 
obeying  N2| N,  M)  =  N(N  +  1)|V,  M)  and  NZ\N,  M)  = 
M\N,  M).  As  E  is  increased,  eigenstates  with  the  same  M 
mix,  forming  superpositions  of  | N,  M)  with  different  N. 
We  denote  the  state  at  finite  E  adiabatically  connected  to 
| N,M)  as  \4>n,m)  anc*  define,  for  notational  conve¬ 
nience,  \N)  =  \(f>N'0),  \N)  =  \(f>N:1),  and  | N)  =  \ 4>N2) 
[see  Fig.  1(a)].  For  appropriate  choices  of  levels,  energy 


(b) 


|i> 


|T> 


FIG.  1  (color  online),  (a)  Eigenenergies  of  H0  =  5N2  —  d(]E 
as  a  function  of  E  (the  top  axis  uses  d  and  B  of  KRb).  Only 
M  >  0  states  are  labeled.  (b),(c)  Examples  of  level  configurations 

{|m0),|mj)}:  (b)  (Va|0)  +  Vl  —  «|I),  |1)>;  (c){Va|2>  +  Vl  —  a\2), 

\fb\  1 )  +  7?|1)  +  Vl  b  —  c|2)}.  Red  (blue)  levels  make  up  the 
dressed  rotor  state  |m0)  (| mj)). 


conservation  ensures  that  dipole-dipole  interactions 
between  two  molecules  preserve  their  total  Nz,  resu¬ 
lting  in  the  many-body  interaction  Hamiltonian  Hdd  = 

3  +  \d+d-  +  d^d^)].  We 

now  focus  on  the  two  configurations  of  levels  and  micro- 
wave  fields  shown  in  Figs.  1(b)  and  1(c),  which  allow, 
respectively,  for  the  realization  of  the  simple  case  Jz  = 

V  =  W  =  0  and  for  the  full  tunability  of  JZ,J±,  V,  and  W. 

In  the  first  configuration,  states  \m0)  =  V«|0)  + 

V I  —  a\\)  and  |nq)  =  |1)  [Fig.  1(b)]  are  chosen,  where 

the  coefficient  a  in  the  dressed  state  \m0)  is  controlled  by 
the  ratio  between  the  Rabi  frequency  and  the  detuning  of  a 
er+  -polarized  microwave  field  acting  on  the  |0)  —  |1)  tran¬ 
sition  [27].  In  40K87Rb,  microwave  intensity  of  a  few 
W / cm2  is  sufficient  to  address  all  hyperfine  levels  equally. 
Projecting  on  states  |0),  |1),  and  |I)  and  keeping  only 
energy-conserving  terms,  we  obtain  d^d^  =  /U-qj  |01)  X 
(10|,-y  +  H.c.  +  and  d(+d^  +  = 

-At^lOlXiOly  -  /ijjlllxnly  +  H.C.,  where  k  G  {i, ./'}, 
5  G  {0,  1, 1},  /uL s  and  /xSJ/  =  (slis^lV)  for  the 

appropriate  p.  The  minus  sign  in  front  of  pL~-{  and  /x 2 (  is 
crucial  to  the  tunability  of  V,  W ,  Jz,  and  J ±  and  appears 
because  two  dipoles  rotating  in  the  x-y  plane  interact  on 
average  attractively  [27].  Projecting  on  |m0)  and  \m , )  and 
again  keeping  only  energy-conserving  terms,  we  obtain 

'^Bp\mpnip)(mpmp\ij  T  ) '  A  p  Aa  \  m  p  ^{nt  p  n2a\  j  j 

p  p,q 

+  :y(|m0mi)<mim0|y  +  H.c.)) 

where  p,  q  G  {0,  1},  A0  =  ap0  +  (1  —  a)pL\,  =  pux, 
Bq  =  — /u,2jfl(  1  -  a),  Bi  =  0,  and  7j_  =  2 - 
/U,2j(l  —  a).  Ap  can  be  thought  of  as  the  effective  dipole 
moment  of  \mp),  while  Bp  is  the  contribution  to  Hdd  from 
transition  dipole  moments  between  rotor  states  making  up 
| mp).  From  the  comparison  of  this  equation  and  Eq.  (1),  we 
can  read  off  V  =  [(A0  +  Ax)2  +  B0  +  B{]/ 4,  IV  =  (A5  + 
B0  —  A2  —  B{)/ 2,  and  Jz  =  (A0  —  Ax)2  +  B0  +  Bx.  A 
simple  case  Jz  =  IV  =  V  =  0  and  J ±  >  0  studied  below 
can  be  approximately  implemented  by  using  a  =  1 ,  which 
does  not  require  a  microwave  field.  A  small  E  field 
(dE/B  >  0. 1  in  40K87Rb)  is  needed  to  prevent  dipole- 
dipole  and  hyperfine  interactions  from  populating 
|01±1).  At  dE/B  =  0.1,  V,  IV,  and  Jz  are  2  orders  of 
magnitude  smaller  than  Jj_  and  can  be  neglected. 

The  second  configuration  we  consider  is  {| m0),  \m , )}  = 
{■sfa\2)  +  Vl  —  a\2),  Vfi|l)  +  V^U)  +  Vl  —  b  —  c|2)} 
[Fig.  1(c)].  The  three  microwave  fields  shown  in  the  figure, 
which  allow  us  to  control  the  coefficients  a,  b,  and  c, 
together  with  the  dc  electric  field  constitute  four  knobs 
that  allow  for  the  full  control  over  the  four  coefficients  V, 
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W,  Jz,  and  J [31].  We  find  that  in  a  small  sphere  in  the 
four-dimensional  ( dE/B ,  a,  b,  c )  space  around  the  point 
(dE/B,  a,  b,  c)  =  (2.97,  0.059,  0.56,  0.38),  where  Jz  = 
7j_  =  V  =  W  =  0,  one  can  achieve  any  value  of  V,  W, 
Jz,  and  J up  to  an  overall  positive  prefactor.  Similarly,  in 
the  special  case  where  a  single  microwave  couples  |I)  and 
1 1),  we  have  a  two-dimensional  (dE/B,  b)  subspace  (with 
a  =  0,  c  =  1  —  b),  in  which  Jz  and  J  ±  can  be  fully 
controlled — for  simulations  of  the  XXZ  model  at  unit 
filling — around  (dE/B,  b)  =  (4.36,  0.56),  where  they 
both  vanish.  While  these  examples  prove  full  controllabil¬ 
ity  in  their  respective  cases,  for  any  desired  relationship 
between  V,  W,  Jz,  and  J±,  there  is  likely  a  different  level 
configuration  that  gives  stronger  interactions  and  uses 
weaker  E,  lower  microwave  intensity,  and/or  more  conve¬ 
nient  microwave  frequencies. 

Phase  diagrams  of  the  nearest-neighbor  and  dipolar 
t-J±  chains. — The  full  tunability  of  the  t-J-V-W  model 
provides  access  to  a  great  variety  of  models  with  poten¬ 
tially  exotic  physics.  As  the  simplest  example  of  this 
physics,  we  present  in  Fig.  2  the  ID  phase  diagram  in  the 
limit  V  =  W  =  Jz  =  0,  which  is  one  of  the  simplest  ex¬ 
perimentally  achievable  cases  (see  above).  Before  analyz¬ 
ing  dipolar  interactions,  we  present  in  Fig.  2(a)  the  phase 
diagram  of  the  nearest-neighbor  t-J ±  chain  as  obtained  by 


FIG.  2  (color  online).  Phase  diagrams  of  the  t-J  ±  chain  with 
(a)  nearest-neighbor  and  (b)  dipolar  J  ]_ .  We  identify  a  metallic 
(repulsive  Luttinger  liquid)  phase  with  dominant  spin-density- 
wave  correlations  (SDW),  a  gapless  superfluid  with  dominant 
triplet  and  singlet  superfluid  correlations  (TS/SS),  a  singlet 
superfluid  with  a  spin  gap  (SS  +  SG),  and  phase  separation 
(PS).  The  spin  gap  is  S  0.35f  in  (a)  and  S  O.lt  in  (b).  Solid 
lines  indicate  phase  transitions  (green,  SG  closes;  red,  inverse 
compressibility  becomes  zero);  dashed  lines  are  extrapolations. 
The  numbers  show  the  value  of  Kp  on  the  dotted  lines.  The  line 
Kp  =  2  is  a  crossover  line  within  the  SS  +  SG  phase  [32].  The 
shaded  region  in  (b)  displays  Kp  =  1  ±  0. 15  as  an  estimate  of 
the  numerical  accuracy. 


using  the  density  matrix  renormalization  group  method  and 
following  the  analysis  of  Ref.  [9].  The  diagram  is  qualita¬ 
tively  similar  to  that  of  the  standard  t-J  chain  [9]:  At 
fillings  n  <  0.65,  we  identify  a  repulsive  Luttinger  liquid 
(Luttinger  parameter  Kp  <  1)  with  dominant  spin-density- 
wave  (SDW)  correlations,  an  attractive  Luttinger  liquid 
(Kp  >  1)  with  dominant  singlet  and  triplet  superfluid  cor¬ 
relations,  a  singlet  superfluid  with  a  spin  gap,  and  phase 
separation.  At  larger  fillings,  the  spin  gap  is  always  zero, 
but  the  other  phases  remain. 

In  Fig.  2(b),  we  present  the  phase  diagram  of  this  system 
in  the  presence  of  dipolar  interactions.  At  low  fillings, 
SDW,  the  gapped  singlet  superfluid,  and  phase  separation 
are  obtained.  This  suggests  that  experiments  with  ultracold 
molecules  can  be  a  useful  tool  for  exploring  the  phase 
diagram  of  the  standard  t-J  model  for  arbitrary  values  of 
J  in  contrast  to  J  <5C  t  in  ultracold  atom  realizations.  At  the 
same  time,  crucially,  both  diagrams  in  Fig.  2  feature  a 
significant  enhancement  of  the  superfluid  region  compared 
to  the  original  t-J  model  because  the  absence  of  attractive 
density-density  interactions  suppresses  phase  separation. 
Furthermore,  the  maximum  value  of  the  spin  gap  in  the 
dipolar  t-J ±  chain  (—0.7 1)  is  twice  that  in  the  nearest- 
neighbor  t-J i  chain  (—  0.35f),  which  is,  in  turn,  twice  that 
in  the  original  t-J  model.  The  larger  spin  gap  should 
facilitate  its  experimental  observation.  Note  that,  as  a  con¬ 
servative  estimate  of  the  numerical  errors  in  computing  Kp 
in  the  presence  of  long-range  interactions  with  up  to  100 
sites,  we  estimate  the  true  line  Kp  =  1  to  lie  in  the  shaded 
region  0.85  <  Kp  <  1.15.  Thus,  the  gapless  superfluid 
cannot  be  identified  in  Fig.  2(b)  within  our  numerical 
precision.  Furthermore,  while  the  line  Kp  =  1  and  the 
line  where  the  spin  gap  closes  may  coincide,  our  analysis 
of  the  correlation  functions  cannot  rule  out  the  existence  of 
an  exotic  intermediate  phase  with  a  spin  gap,  Kp  <  1 ,  and 
dominant  superfluid  correlations. 

Preparation  and  detection. — Ground  states  at  specific 
points  in  the  phase  diagram  can  be  prepared,  e.g.,  by 
applying  an  additional  microwave  field  coupling  |m0) 
and  | mx)  and  performing  an  adiabatic  passage  from  an 
easily  accessible  state  to  the  desired  ground  state  by  tuning 
the  Rabi  frequency  and  the  detuning  of  the  microwave  field 
[18].  Direct  probing  of  molecules  [33]  and  conversion  of 
molecules  back  to  atoms  [5]  can  in  principle  both  be  used 
for  detection  via  noise  correlations  in  the  time  of  flight  [1] 
or  via  in  situ  single-site  imaging  [34,35]. 

As  a  specific  example  of  a  detection  technique  available 
in  current  experiments,  we  propose  to  use  Bloch  oscilla¬ 
tions  [36,37].  In  Fig.  3,  we  present  our  results  as  obtained 
via  the  Krylov-space  variant  of  the  adaptive  time- 
dependent  density  matrix  renormalization  group  method 
[38,39]  for  20  sites  when  adding  a  linear  field  along  the 
chain  at  filling  n  =  0. 1  and  n  =  0.2  for  different  values  of 
J i_/t.  In  the  singlet  superfluid,  the  amplitude  of  the  oscil¬ 
lations  drops  and  the  frequency  of  the  oscillations  doubles 
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FIG.  3  (color  online).  Bloch  oscillations  in  a  dipolar  t-J ± 
chain  of  20  sites  with  open  boundary  conditions  at  fillings 
n  =  0. 1  and  0.2,  and  tilting  field  E&J t  =  1  per  site.  The  main 
plot  shows  the  difference  between  the  center-of-mass  position 
jcCm  at  time  zero  and  at  the  first  minimum.  The  vertical  line 
indicates  the  approximate  value  of  J  ±  / t,  at  which  the  spin  gap 
closes  for  n  =  0.2.  The  inset  shows  the  time  evolution  of  xCM  at 
n  =  0.2  at  three  indicated  values  of  J±/t. 

(see  7_l /t  =  7  in  the  inset)  relative  to  the  gapless  phase  due 
to  the  presence  of  bound  pairs.  Bloch  oscillations  should  be 
observable  in  direct  absorption  imaging  [33].  While  neither 
the  frequency  nor  the  amplitude  of  the  oscillations  show 
any  sharp  features  at  this  small  system  size,  the  fit  of 
experimental  data  to  numerical  results  should  allow  for 
the  location  of  the  phase  transition  (vertical  line  in  Fig.  3) 
even  for  small  system  sizes.  As  a  complementary  method 
for  identifying  the  transition,  we  propose  spectroscopic 
measurement  of  the  spin  gap  [40]. 

Outlook. — We  have  presented  a  toolbox  for  simu¬ 
lating  a  highly  tunable  anisotropic  t-J-V-W  model  with 
polar  molecules.  The  advantages  of  this  molecular  toolbox 
over  its  atomic  counterpart  are  higher  energy  scales  and 
independent  tunability  of  interactions  and  tunneling.  This 
toolbox  should  enable  the  simulation  of  condensed  matter 
phenomena,  as  well  as  the  stabilization  and  controllable 
preparation  of  unconventional  phases,  such  as  d- wave 
superfluids  [3].  The  phase  diagram  of  the  experimentally 
simple  case  of  a  t-J±  chain  shows  an  enhanced 
superfluid  region,  which  we  propose  to  probe  via  Bloch 
oscillations.  We  expect  that  the  2D  t-J-V-W  model 
can  similarly  be  tuned  into  exhibiting  enhanced 
superfluidity  [41]. 

The  present  Letter  also  opens  other  exciting  research 
avenues.  In  particular,  natural  extensions  of  the  model 
include  [30]  spatially  anisotropic  interactions  produced 
by  a  tilt  in  the  dc  electric  field,  spin-dependent  tunneling 
obtained  by  adjusting  lattice  beams,  S  >  1/2  models  real¬ 
ized  by  choosing  more  than  two  dressed  states,  and  systems 
with  an  orbital  degree  of  freedom  encoded  in  the  nuclear 
spin.  Moreover,  by  considering  molecular  Wigner  crystals 
[25,42,43],  where  the  intermolecular  distances  are 
smaller  than  in  an  optical  lattice,  one  can  further  increase 
the  interaction  strength.  Furthermore,  by  analogy  with 
Ref.  [18],  we  expect  our  ideas  to  be  extendable  to 


Rydberg  atoms.  Finally,  one  can  envision  applications  of 
the  present  system  to  quantum  computation  (especially 
if  one  uses  nuclear  spin  to  store  information),  precision 
measurements,  and  controlled  quantum  chemistry  [4]. 
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